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Abstract 

Calibration of polarimetric imaging SARs using point calibration targets is discussed. 
The four-port network calibration technique [5] is used to describe the radar error model. 
The polarimetric ambiguity function of the SAR is then found using a single point tar- 
get, namely a trihedral corner reflector. Based on this, an estimate for the backscattering 
coefficient of the terrain is found by a deconvolution process. 

A radar image taken by the JPL aircraft SAR is used for verification of the deconvolution 
calibration method. The calibrated responses of point targets in the image are compared 
both with theory and the POLCAL technique [8]. Also, responses of a distributed target 
are compared using the deconvolution and POLCAL techniques. 

1 Introduction 

Calibration techniques available in the literature can be categorized into two major 
groups: 1) calibration techniques for imaging radars, and 2) calibration techniques for point- 
target measurement systems, which may also be appropriate for imaging radars. In the first 
group, the scattering properties of clutter are usually employed to simplify the calibration 
problem such as in Sheen and Kasischke [1989], van Zyl [1990], and Klein [1989]. Among the 
point-target calibration techniques, are Whitt, et al. [1990], Barnes [1986], and Sarabandi, 
et al. [1990] which uses a sphere and any other depolarizing calibration target. However, the 
isolated antenna assumption can lead to significant errors when the ratio of cross- to like- 
polarized terms is small and/or cross-talk contamination is large. To remove this drawback, 
the single target calibration technique (STCT) has recently been developed [Sarabandi and 
Ulaby, 1990]. 

The main thrust of this paper is to show how the point target calibration techniques can 
be applied to imaging SARs. In particular, the STCT will be employed here since it only 
requires one calibration target. This method is then compared with the POLCAL technique 
[8] which has been developed specifically for imaging radars. 

2 A New Approach for Calibration of Imaging SARs 

In this new method the polarimetric ambiguity function of the SAR processor as well 
as the distortion matrices of the radar system are found from a trihedral corner reflector 
response. First a summary of the single target calibration technique (STCT) is given and 
then a theoretical model which relates the point target response to distributed targets will be 
developed. The model in conjunction with the STCT is then used to obtain a deconvolution 
matrix for estimation of the backscattering coefficient. 
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2.1 Single Target Calibration Technique 

In this technique the antenna system and two orthogonal directions in free space are 
modeled as a four-port passive device. The measured scattering matrix of a target with real 
scattering matrix s is approximated by 
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The complex quantities /? p , T q are the receive and transmit channel distortions with p,q = 
u, h and C is the antenna cross-talk factor. Once the distortion parameters of the radar 
system are found through the measurement of a trihedral corner reflector, U t , with radar 
cross section <J t whose scattering matrix is diagonal, the actual scattering matrix of the 
target can be obtained from 
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The single target calibration technique (STCT) has been tested both under laboratory and 
field conditions using L-, C-, and X-band scatterometer systems and it has been shown that 
a calibration accuracy of 0.5 dB in amplitude and 5 degrees in phase can be achieved. 


2.2 Difficulties in Calibration of SARs Using Point Targets 

To understand the steps involved in calibration of imaging SARs using point calibration 
targets, the generation of a high resolution image from raw data must be examined. The 
received raw data in each of the channels of a polarimetric SAR can be described by 

Upq (<) = J A Spg i x ', y') f (< - dx'dy' p,q = v,h (6) 

where A is the illumination area by the physical antenna, ( x ', y ') is the reflectivity of 
the terrain being mapped (cTp q (x,y) = 47r | Sp q (x,y) | 2 ), and p and q are the polarization 
state of the receiver and transmitter respectively. Function / ( t ) is a particular wave form 
radiated by the transmitter and can usually be represented by 

f{t)=g{t)e^ (7) 
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where g(t ) is a slowly varying function and u> 0 is the angular frequency of the radar system. 
In equation (6) R is the distance from the antenna to the scattering point (x', y') on the 
ground. 

One way to retrieve the backscattering coefficient o® q (x,y) from the received signal 
U pq ( t ) is to pass the signal through a matched filter having an impulse response /* (t — 
[Cutiona, 1970], that is 

U 'pq ( x . V) = JJ A a P 9 ( x '> y') / ^ /* dx'dy'dt (8) 

By performing the integration with respect to time, the quantity 

(*, y; y ') = Jj^t- f* (j - dt (9) 

known as the ambiguity function can be obtained. Therefore (8) can be written as 

U' pq {x,y)= \ a 0 pq (x,y')xl>(x,y,x',y')dx'dy' . (10) 

J A 

If the ambiguity function is a Dirac delta function, i.e. 0(x, y; x', y') — 8{x — x',y — ?/'), the 
backscattering coefficient can be directly obtained from (10). 

By substituting (7) into (9) and assuming that g(t) is a linearly frequency-modulated 

/ . 7tA/ ^2\ 

pulse of duration r g(t ) = e t and that the integration time in (9) is over N + 1 
pulses of the transmitter we have 
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Using (11) in (10) does not resolve cr® q (x,y) completely. To obtain a better estimate a 
deconvolution process must be attempted. This requires that the calibration technique be 
able to estimate the ambiguity function. 


2.3 Calibration Procedure and Estimation of Backscattering 
Coefficient 


The error model for a polarimetric SAR must include the uncertainties such as the 
antenna cross talk and channel imbalances as well as the uncertainties in the ambiguity 
function. Suppose the radar system is linear. The polarimetric response of the terrain with 
polarimetric reflectivity S° (x,y) is given by 


U(x,y) = R [ S° (x 7 , y') ip(x, y; x', y')dx'dy' 
.J A 


( 12 ) 


In (12) the amplitude and phase of the propagation factor ^ ^ ^ has been excluded 
and R and T are, respectively, the receive and transmit distortion matrices. Assume 
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xp (x, t/; x ', y ') = ip (x — x', y — y'). Since radar images are discretized into a finite number of 
pixels, the discretized form of (12) must be considered, where the integral is approximated 
by a double summation, thus 


U (m, n) = R 


S°(i,j)ip{m-i,n-j)AxAy 
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By changing the index of the summations such that the maximum of the ambiguity func- 
tion occurs at (0,0), and lumping together the ambiguity function and distortion matrices, 
equation (13) becomes 


N_ _M 

U(m,n) = AxAy ^ ^ AR :J S° (m — z, n — j) AT (14) 



where AR tJ and AT q are the receive and transmit ambiguity-distortion matrices respec- 
tively. 

If a trihedral with radar cross section crj located at (xo, yo) is used as a calibration target, 
its reflectivity function is expressed by 


Sj (x, y) 
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where I is a 2 x 2 identity matrix. Substituting (15) into (12) results in the polarimetric 
response of the SAR to the trihedral. When the discretized form of this response is combined 
with the resultant ambiguity function and distortion matrices, the response becomes 
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where i = m — k and j = n — l for a trihedral at the A:/th pixel. Figure (1) shows the measured 
amplitude of the polarimetric response of the JPL aircraft SAR to an 8-foot trihedral. 

An approximate form of (14) can be obtained by assuming that over the pixels where 
the ambiguity function is non-zero the reflectivity of the terrain is constant and is equal to 
that of the center pixel, giving 
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If the four-vector form of U (m, n) and S° (m, n) are, respectively, denoted by 
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then (17) becomes 
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where D is the deconvolution matrix which is independent of pixel coordinates m and n. 
Therefore the calibrated estimate of the reflectivity matrix is given by 


- o 


«S (m,n) = D 1 U (m,n) . 


3 JPL POLCAL Technique 

POLCAL [van Zyl et al., 1990] is a combination of phase calibration by Zebker and 
Lou [1990] and a calibration algorithm by van Zyl [1990] based on properties of distributed 
targets. In this technique the radar error model and corrections are done in three steps. 
The first step is phase calibration where the radar distortion matrices are assumed to be 
diagonal with only phase differences. That is, the measured scattering matrix is assumed 
to be 
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where s pq are the theoretical values for the scattering matrix elements. Reciprocity mandates 
that Shv = Svh and therefore the quantity Uh v ul h must have zero phase. The phase difference 
< fa — <f> r is averaged over the entire image, then subtracted from u / lv , to form a matrix Z 
whose off-diagonal elements have almost identical phases. Next this matrix is symmetrized 
by averaging the off-diagonal elements and then stored in the matrix Y. The data is then 
coded and stored in the form of the Stokes scattering operator, with groups of four adjacent 
pixels in a row being summed. The quantity <f> T + <f> t is obtained from a trihedral response 
by calculating Y hh Y* v . 

The next steps include cross-talk removal and adjusting for co-channel gain balance. Here 
the radar error model of the phase calibrated symmetrized response (Y') is represented by 
reciprocal transmit and receive distortion matrices, i.e. 
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where s is the actual scattering matrix of the target, and 62 represent the antenna cross- 
talk, and / is the co-channel imbalance. Using the above equations <$1 and 62 / f are found 
iteratively. The amplitude of the co-channel imbalance / is obtained by calculating the ratio 
of the total power of VV to HH of a corner reflector response over 16 x 16 surrounding pixels. 
However, it is left unjustified why the HH and VV responses should be added noncoherently. 


4 Results and Comparison 

This section presents the results of applying the three different calibration techniques 
described in sections 2.1, 2.3, and 3 to the same scene. Each technique was applied twice, 


171 



each time with a different calibration target, which helped in determining the reproducibility 
of the calibration. 

4.1 Data Formats and Test Scene 

The results given in this section were obtained by processing the same JPL AirSAR scene 
as provided by JPL in two different formats. The format used by the techniques described in 
section 2 is the so-called “hires” format, which provides the four scattering matrix elements 
as single-precision (4-byte) complex numbers. The format used by the POLCAL technique is 
radically different. The symmetrization, quantization from 4 bytes to 1, and summing from 
1-look to 4-look, as described in section 3, have been carried out to produce the so-called 
“compressed” format. 

Figure 2 shows the total power image of the particular scene used during this study. 
Figure 3 shows the uncalibrated HH polarized raw data in the vicinity of the three trihedrals. 
Note that the trihedrals are identical but that the responses are not, which is of great concern 
when calibration depends on identical responses to identical targets. Also, a distributed 
target, as outlined in Fig. 2, was used for comparisons despite the lack of any known values 
for its cross sections. 

4.2 Assessment of STCT 

Figure 4 shows the calibrated polarization signatures for the trihedral that was not used 
as calibration target, in each of the two cases. The results are excellent: the levels are within 
1 dB of theory, with the co-channel imbalance no more than 0.5 dB. Also, the cross-to-like 
isolation has improved from 20 dB before calibration to about 40 dB. 

Figure 5 shows the calibrated signatures of two different PARCs: VV and 45°. The VV 
PARC signature is as expected with a 23 dB isolation between VV and HH. The 45° PARC 
signature has the peak very nearly centered on x = 0, $ = 45°, as expected. Because of 
this favorable comparison with the expected results one can conclude that this technique is 
applicable to the JPL AirSAR imaging radar data. 

4.3 Assesment of POLCAL 

Because the calibration trihedral is used only to determine the co-channel imbalance and 
absolute level, this calibration technique gives noticeably different results when compared 
to those in section 4.2. Figure 6 shows the single-pixel signatures of various trihedrals after 
calibration with POLCAL. Here the co-channel imbalance is between 0.5 dB and 1 dB, 
slightly worse than for STCT. The cross-to-like isolation is generally 23 dB, 17 dB worse 
than for STCT, while one has an isolation of 175 dB, apparently a fluke since it appears 
only once. 

To compare the absolute levels, the signatures are computed by summing over the same 
region that is used to sum the powers of the calibration target. These results are shown in 
Fig. 7. Generally the signatures do not as closely resemble a trihedral as did those when 
using just a single pixel (Fig. 6). The absolute levels, however, are within 1 dB of the 
expected levels. 

Figure 8 shows the same two PARCs as before, however they look quite different using 
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the POLCAL calibration scheme. The VV PARC has a 15 dB isolation, 8 dB worse than 
STCT, while the cross-pol signature is significantly distorted due to the symmetrization 
step: the 45° PARC now has two cross-pol peaks, at ±45°. 

4.4 Comparison using a distributed target 

First, a comparison of the two formats, uncalibrated, is in order. Figure 9 shows that 
the hires data has a VV power of about 3.75 dB less than HH, while the compressed data 
has them 4.75 dB apart. There are similar differences after calibration as well. 

The deconvolution calibration scheme described in Sec. 2.3 was used, with the summation 
on D,j done over all pixels in the vicinity of the calibration target that had a cross-pol 
magnitude about 20 dB above the noise. The results are shown in Fig. 10. The signatures 
are very similar, the only difference being in their absolute levels - a difference of at most 0.75 
dB. Hence, the differences in the measured ambiguity function for the different calibration 
targets do not significantly affect the calibration of distributed targets. Figure 11 shows the 
results of the POLCAL calibration and looks very similar to the results in Fig. 10. The 
difference between the two calibration schemes is only on the order of 1 dB. 

5 Conclusions 

A new method for calibration of polarimetric imaging SARs using point targets has been 
developed. The technique requires a single calibration target, namely a trihedral, to find 
the ambiguity-distortion matrix (polarimetric ambiguity) function of the SAR system. 

The validity of the technique is examined by calibrating a scene which includes a variety 
of point targets with known scattering matrices. It is found that the error model provided 
by STCT is an appropriate one and enhances the measured polarization of the point targets. 
The deconvolution technique compares favorably with the POLCAL technique for absolute 
radiometric calibration of distributed targets. This new technique is preferred to POLCAL 
only if the radar system has large distortions. 
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Figure 1: Measured Ambiguity-Distortion 
Matrix values in dB for a single trihedral. A 
small region surrounding the brightest pixel 
is shown. 

Figure 3: Measured HH-polarized values 

in dB for all three trihedrals. A small re- 
gion surrounding them is shown. Note that 
the response is significantly different for the 
identical trihedrals. 



Figure 2: Total Power image of calibra- 

tion targets and surrounding area. Shown 
are all three trihedrals, the five PARCs, and 
the distributed target area to the left. 
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Power values: hh* 32.72 dB. vv* 32.19 dB, bv*-14.77 dB, vb*-12.12dB 
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Power values: hh* 34.66 dB, vv* 35.19 dB. bv- -11.81 dB. vh= -8.97dB 



Figure 4: Validation results for STCT. (a) 
shows trihedral #2 when the calibration tar- 
get was trihedral #1. (b) shows trihedral 

#1 when the calibration target was trihedral 
# 2 . 


PoaeroikBK bfc* -5 07 dB. vv* 17.(7 dB Po»er »alu»: hb= 16 II dB. vv* 14 61 dB 
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bv- -6.90 dB. vb*-29 76 dB bv- 1731 dB. vh* 17 10 dB 



Co Pol Crow Pd Co Pol Crow Pol 


VV PARC (b) 45° PARC 

Figure 5: Validation results for STCT. (a) 
shows the VV and 45° PARCs when the cal- 
ibration target was trihedral #1. (b) shows 
the VV and 45° PARCs when the calibration 
target was trihedral #2. 



Power value* t*= 2130 dB. w- 28.64 dB 
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bv* 5.85 dB 


trihedral #1 (a) trihedral #2 
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Power value* to* 2839 dB. w- 28.94 dB 
bv- 4.64 dB 


fab- 29 43 dB, w- 3033 dB 
bv- 14337 dB ^ 


trihedral #1 (b) trihedral #2 trihedral # 1 (b) trihedral #2 

Figure 6: Validation results for POLCAL, Figure 7: Validation results for POLCA 
displaying only 1 pixel, (a) shows both tri- displaying sum over adjacent area, (a) sho 1 
hedrals when the calibration target was tri- both trihedrals when the calibration targ 
hedral #1. (b) shows both trihedrals when was trihedral #1. (b) shows both trihedr; 
the calibration target was trihedral #2. when the calibration target was trihedral $ 
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Power value* bb = 6 *0 dB, w= 21.64 dB Power values bh = 21.96 dB. w* 17.72 dB 
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VV PARC (b) 45° PARC 

Figure 8: Validation results for POLCAL, 
displaying only 1 pixel, (a) shows the VV 
and 45° PARCs when the calibration target 
was trihedral #1. (b) shows the VV and 

45° PARCs when the calibration target was 
trihedral #2. 
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Figure 9: Uncalibrated Signature of dis- 

tributed target region, (a) shows the hires 
format, (b) shows the compressed format. 


Power values: hh* -7.12 dB, vv* -10.95 dB, bv= -14.53 dB, vb* -15.17 dB 
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P ower values: hh* -6.99 dB, vv* -10.60 dB, bv* -14.52 dB, vb* -14.71 dB 



Figure 10: Signature of distributed target 
region, Deconvolution calibration, (a) shows 
the response with trihedral #1 as the cali- 
bration target, (a) shows the response with 
trihedral #2 as the calibration target. 


Power values: hb* -7.19 dB, vv* -1 1.29 dB, bv* -15.06 dB 



Power values: hb* -7.31 dB, vv*-10.9«dB. bv* -14.96 dB 



Figure 11: Signature of distributed target 
region, POLCAL calibration, (a) shows the 
response with trihedral #1 as the calibration 
target, (a) shows the response with trihedral 
#2 as the calibration target. 
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